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We develop an atomistic, nearest-neighbor sp 3 s* tight-binding Hamiltonian to investigate the 
electronic structure of dilute bismide alloys of GaP and GaAs. Using this model we calculate that 
the incorporation of dilute concentrations of Bi in GaP introduces Bi-related defect states in the 
band gap, which interact with the host matrix valence band edge via a Bi composition dependent 
band anti-crossing (BAC) interaction. By extending this analysis to GaBi x Asi-a, we demonstrate 
that the observed strong variation of the band gap (E g ) and spin-orbit-splitting energy (A so) 
with Bi composition can be well explained in terms of a BAC interaction between the extended 
states of the GaAs valence band edge and highly localized Bi-related defect states lying in the 
valence band, with the change in E g also having a significant contribution from a conventional 
alloy reduction in the conduction band edge energy. Our calculated values of E g and Aso are 
in good agreement with experiment throughout the investigated composition range (x < 13%). 
In particular, our calculations reproduce the experimentally observed crossover to an E g < Aso 
regime at approximately 10.5% Bi composition in bulk GaBi^Asi-^. Recent x-ray spectroscopy 
measurements have indicated the presence of Bi pairs and clusters even for Bi compositions as low 
as 2%. We include a systematic study of different Bi nearest-neighbor environments in the alloy to 
achieve a quantitative understanding of the effect of Bi pairing and clustering on the GaBi^Asi-x 
electronic structure. 

PACS numbers: 61.43.Dq, 71.15.-m, 71.20.Nr, 71.22.+i, 71.23.-k, 71.23.An, 71.55.-i, 71.55.Eq 



I. INTRODUCTION 

Highly mismatched semiconductor alloys such as 
GaN^Asi-z have attracted considerable interest both 
from a fundamental perspective and also because of their 
potential device applications^^. When a small fraction 
of As is replaced by N in GaAs, the band gap (E g ) ini- 
tially decreases rapidly, by sa 150 meV when 1% of As 
is replaced by INPEl. Similar behavior has been experi- 
mentally observed in GaBi^Asi-a;, where EL has been 
observed to initially decrease by 90 meWES! per % of 
Bi replacing As in th e alloy. In addition to this, re- 
cent experiment d 11 * 12 ! have also revealed the presence of 
a large bowing of the spin-orbit-splitting (SO) energy 
(Ago) with increasing Bi composition. 

It has been recently showrP^ that by increasing the 
Bi composition in GaBi x Asi_ x to « 11% we enter an 
Eg < Ago regime in the alloy. This regime is of interest 
for the design of highly efficient optoelectronic devices 
since it opens up the possibility of suppression of the 
non-radiative CHSH Auger recombination process, a loss 
mechanism which plagues the efficiency and dominates 
the threshold current in III-V multinary lasers operating 
in the telecommunication wavelength ranged 

The band gap reduction in GaN-rAsi-z has b een ex- 
plained using a band anti-crossing (BAC) models It 
is well established that replacing a single As atom by N 
introduces a resonant defect level above the conduction 
band edge (CBE) in GaAs. This occurs because N is con- 
siderably more electronegative and is also 40% smaller 
than As. The interaction between the resonant N states 
and the CBE of the host GaAs matrix then accounts for 



the observed rapid reduction in E g . 

It has been proposed that replacing As by Bi could 
introduce a similar BAC effeclP. Like N and As, Bi is 
also a group V element, which lies below Sb in the pe- 
riodic table. Bi is much larger than As, and is also less 
electronegative. It should therefore be expected that any 
Bi-related resonant defect levels should lie in the valence 
band (VB) and that, if an anti-crossing interaction oc- 
curs, it will occur between the Bi-related defect level and 
the valence band edge (VBE) of the GaAs matrix. 

Recent x-ray absorption spectroscopy analysis of 
GaBia;Asi_J^ revealed that Bi atoms tend to form pairs 
and clusters in the alloy for x > 2%. For this reason it 
is also of interest to explore the impact of Bi pairs and 
Bi clusters on the electronic structure of GaBi^Asi^. 
Such analysis requires atomistic modelling of the alloy 
since continuum methods such as the k • p and effective 
mass approximations are unable to predict the influence 
of such pairs and clusters on the electronic structure. 

Although the re is an increasing e xperime ntal interest 
in the structura l 14 * 15 ! and electronic ^ 9 ! 11 ! 16 ] properties of 
GaBi^Asi-a;, there have been comparatively few theo- 
retical studies to date of this extreme semiconductor al- 
loy. Previo us th eoretical investigations have employed 
the LDA+d 17 * 18 ! , k • pP, and empirical pseudopotentiaP^ 
methods. These theoretical approaches have provided 
a range of insights, but it remains controversial as to 
whether the observed reduction in energy gap is due to 
an anti-crossing interaction^ or is better explained using 
more conventional alloy and disorder models^. 

The LDA+C calculations of Janotti et alP^ presented 
the first band structure calculations of GaBi, and ini- 
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FIG. 1: Bulk band structure of GaBi. The 
broken red lines are the bands obtained 
from our sp 3 s* parametrisation. The solid 
black lines are the LDA+C bands calcu- 
lated by the authors of Ref. 11. (used with 
permission) 



tial investigations of the influence of Bi on the electronic 
structure of GaBi^Asi^a;. Zhang et alW^ presented pseu- 
dopotential calculations which were in good agreement 
with the experimentally measured values of E g and A50 
over the small composition range considered (x < 4%). 
They also identified a Bi-related resonant defect state, 
Eb%, which they found to lie 80 meV below the GaAs 
VBE in the dilute doping limit. The best choice of model 
to describe the influence of such a resonant level on the 
band structure has recently been discussed by Deng et 
aZP^, who argued that it was not appropriate to use the 
BAC model to describe the impact of Bi incorporation, 
and proposed instead the use of a band-broadening pic- 
ture. 

We find here a Bi-related level below the GaAs VBE 
in GaBij-Asi-a;, in agreement with previous work. Our 
comparative calculations of GaBi^Pi-a; and GaBi^Asi-^ 
confirm an anti-crossing interaction between the Bi- 
related resonant state and the VBE. We also investi- 
gate the influence of Bi pair and cluster states on the 
electronic structure, showing that these pairs and clus- 
ters give rise to additional energy levels which tend to 
move into the energy gap with increasing cluster size, 
analogous to the behavior of N pair and cluster levels in 
GaAsPand GaP^. 

The k • p calculations of Alberi et aZP include explic- 
itly BAC interactions between the GaAs VBE and lower- 
lying Bi-related levels in order to fit the experimentally 
observed variation of E g for dilute Bi compositions. They 
also propose a second BAC interaction between the GaAs 
SO band and a Bi-related SO defect level, 1.9 eV below 
the GaAs VBE. The inclusion of such a level would nec- 
essarily produce a BAC-induced upward shift in the SO 
band edge energy, Eso- We find in our calculations that 
the evolution of the spin-split-off band is well described 
using a conventional alloy model, without the need to 
include such an SO-related BAC interaction. 

We have previously demonstrated that a parametrised 



tight-binding (TB) approach provides an accurate and 
versatile method to analyse the electronic structure 
of highly mismatched semiconductor alloys such as 
GaN^Asi-a;, GaBi^Pi-a;, and GaBi^Asi-a;. Since a ba- 
sis of localized states is employed, it is possible using the 
TB method to probe explicitly and quantitatively the ef- 
fect of an isolated impurity atom (e.g. replacing a single 
As atom by either N or Bi in GaAs) on the host ma- 
trix electronic structure^EH. The sp 3 s* TB Hamiltonian 
developed for dilute nitride alloys 2 ^ fully accounts for a 
wide range of experimental data, and also gives results 
in good agreement with pseudopotential calculations^. 
We present here an sp 3 s* TB Hamiltonian to calculate 
the electronic structure of dilute bismide alloys of GaP 
and GaAs. This Hamiltonian gives good agreement with 
pseudopotential and LDA+C calculations 1 -, and with 
experiment over the full composition range investigated 
(x < 13%). 

We start our calculations with the study of ordered 
GaBia-As^a; and GaBi 3 ,P 1 _ :E supercells. It is useful to 
include GaBi^Pi_ x in our analysis because Bi is suffi- 
ciently different to P that it introduces defect states just 
above the VBE in GaP^l. Since these Bi-related defect 
states are in the energy gap, it is then much easier to 
identify and follow the interaction between the VBE and 
these levels as a function of increasing Bi composition. 
Our calculations show clearly for GaBi^Pi-a; that there is 
an anti-crossing interaction between the Bi defect states 
and the GaP VBhP. 

Having established that there is a BAC interaction in 
GaBi^Pi^a;, we then extend our calculations to consider 
ordered supercells of GaBi^Asi-^. By comparison with 
GaBi^Px^a;, we identify that there are Bi-related reso- 
nant defect levels below the VBE in GaBi^Asi-a; and 
that these resonant defect levels interact with and push 
the VBE up in energy through an anti-crossing interac- 
tion. The effect of the resonant levels has been confirmed 
in GaN x Asi_:i; through photoreflectance (PR) measure- 
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merits, which show a higher energy feature . E+ , asso- 
ciated with transitions to the resonant leve l 5 * 26 * ^ . We 
show that, because of the large density of valence states, 
it is not possible to experimentally identify similar higher- 
energy transitions at any Bi composition in GaBLjAsi-a,. 

Next, we investigate the effect of Bi pairs and clusters 
in GaBi-rAsi-j; alloys, by placing a single pair or cluster 
in a large (4096 atom) supercell in order to investigate 
its effect on the electronic structure. Our calculations 
highlight the evolution of the defect energy levels and 
their interaction with the VBE with increasing cluster 
size. This analysis later aids our understanding of the 
calculated electronic properties of randomly disordered 
GaBi^Asi-a, supercells, in which the number of Bi pairs 
increases as x 2 and the number of rt-atom clusters in- 
creases as x n . 

Having established an understanding of the evolution 
of the Bi resonant state in ordered supercells contain- 
ing only a single Bi atom, single Bi pair, or small Bi 
cluster, we next present calculations of disordered 4096 
atom GaBij, Asi^-e supercells with larger Bi compositions 
(x < 13%). In such systems, we find that the CBE en- 
ergy decreases both linearly and rapidly with increasing 
Bi composition, so that the experimentally observed re- 
duction in the band gap is then due both to a BAC in- 
teraction in the VB and to conventional alloy behavior 
in the conduction band (CB). In addition, our calcula- 
tions indicate tha t the observed strong bowing of the 
SO-splitting energjlilMI j s predominantly attributable to 
the upward shift in the GaBi^Asi-^ VBE with Bi compo- 
sition, with a minor contribution from a downward shift 
in Ego, which is linear in Bi composition and hence well 
understood as a conventional alloying effect. This linear 
variation in Eso is contrary to the SO-BAC interaction 
included in the k • p Hamiltonian introduced by Alberi 
et a/P and shows that the Bi-related resonant state in- 
troduced in the VB interacts solely with the GaAs VBE 
at the Brillouin zone center (r point). 

Our calculations accurately reproduce the experimen- 
tally observed^ crossing of the band gap and the SO- 
splitting energies at a Bi composition of w 10.5% in the 
alloy. Our analysis therefore lends support to the experi- 
mental conclusion that GaBi^Asi-^, alloys can be used to 
obtain emission in the technologically important 1.3/xm 
and 1.5/^m wavelength ranges on a GaAs substrate, while 
also demonstrating that the unusual material behavior 
present in the alloy can be understood simply in terms 
of a BAC interaction. 

The remainder of this paper is organised as follows: 
In section II we give a detailed presentation of the TB 
Hamiltonian and calculation procedure. In section III-A 
we apply the TB model to explore in detail the elec- 
tronic structure of dilute ordered bismide alloys of GaP 
and GaAs, beginning with GaBi^Pi^a,. Following this, 
we continue our analysis of GaBi-rAsi-^ by performing a 
systematic study of the effects of Bi pairing and cluster- 
ing on the electronic structure in section III-B. We then 
consider disordered GaBi^Pi^a, supercells in section III- 



C and disordered GaBi^Asi-^ supercells in section III-D. 
Finally we present our conclusions in section IV. 

II. METHODOLOGY 
A. sp 3 s* tight-binding model 

We use an sp 3 s* nearest-neighbor tight-binding Hamil- 
tonian including spin-orbit coupling to investigate the 
electronic structure of the different alloys. The inter- 
atomic interaction parameters vary with bond length and 
the magnitude of the on-site parameters depends explic- 
itly on the overall neighbor environment!^. 

The alloy band structure calculations are under- 
taken using the NanoElectronic Modelling (NEMO 3- 
D) simulatoi | 22HSil. NEMO 3-D is a fully atomistic simu- 
lator and its capabilities include modelling of realistic 
multi-million atom geometries for embedded quantum 
dot devices^H23 ; strained disordered quantum wells^, 
and single impurities in Si FinFET structures^. 

In our sp 3 s* model, the on-site orbital energies for a 
given atom are determined by taking the average of the 
values from the binary compounds formed by that atom 
and its four nearest neighbors, taking into account the 
relevant valence band offsets (VBOs). The inter-atomic 
interaction energies are taken to vary with bond length, 
a, as (^)", where d is the equilibrium bond length in 
the equivalent binary compound, and 77 is a scaling pa- 
rameter whose magnitude depends on the type of inter- 
action being considered. The scaling indices for GaP and 
GaAs are determined by fitting the bond length depen- 
dence of their nearest-neighbor interaction parameters to 
a range of hydrostatic deformation potentials^. Since 
there are no experimental or theoretical deformation po- 
tentials available for GaBi, the scaling indices for GaBi 
were adapted from those of GaAs. We modify the GaAs 
r\ values, in particular for the p states, in order to ensure 
that Bi gives a defect state in the energy gap of GaP, 
as observed experimentalljPH. These modified values of 
77 for GaBi then give a Bi resonant state inside the va- 
lence band of GaAs for GaBiAs, consistent with earlier 
pseudopotential calculations^. The effects of changes in 
bond angle are taken into account through the two-center 
integral expressions of Slater and Kostei^. 

The self-energies and inter-atomic interaction parame- 
ters for the different binary compounds are listed in Ta- 
ble |T] The parameters for GaAs and GaP are based 
on those which we have used previously when study- 
ing GaN^Asi-J^l and GaN^Pi-J^, but we now include 
spin-orbit interactions. The parameters for GaBi were 
obtained by fitting to LDA+C calculations^. Fig. [l] 
shows the bulk GaBi band structure obtained from our 
sp 3 s* Hamiltonian, compared with the LDA+C calcula- 
tions of Janotti et af^l. 

The valence band offsets for GaAs and GaP are taken 
from Ref. 38. Since the binary compound GaBi has yet 
to be isolated, we need to estimate a value for the GaBi 
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FIG. 2: Valence band offset (VBO) as a function of lattice pa- 
rameter for Ga-containing III-V semiconductor binaries. The 
values for GaN, GaP, GaAs and GaSb are taken from Ref. 
38. The GaBi VBO is estimated by extrapolating the trend 
to the GaBi lattice constant (6.328 Af^. For comparison, the 
GaBi VBO value of 0.8 eV estimated by the authors of Ref. 
7 is indicated by a black dot. 



VBO. A VBO of l.leV for GaBi was estimated by an ex- 
trapolation of the VBO as a function of lattice parameter 
(a) for Ga-containing III-V semiconductor binaries (see 
Fig. [5]). This differs from the value of 0.8eV used earlier 
by Alberi et alQ which appears to have been estimated 
by a linear extrapolation of VBO against the lattice pa- 
rameter for GaN and GaP only. The scaling exponents, 
77, used to describe the variation of the inter-atomic ma- 
trix elements with bond length are listed in Table [Tl] 

The TB calculations were carried out on supercells con- 
taining up to 4096 atoms using the NEMO 3-D simula- 
tor. For the disordered structures, Ga204sBij^Y2048-M 
(Y = P, As) supercells were constructed by substituting 
Bi atoms at random Y sites on a GaY lattice and relax- 
ing the crystal to its lowest energy configuration using 
a parametrised valence force field (VFF) model includ- 
ing anharmonic corrections^ to the Keating potential. 
The VFF parameters used are listed in Table |III| The pa- 
rameters for GaAs are those obtained by LazarenkovaS^l, 
while those for GaP and GaBi were derived based on the 
expressions of MartirP2l and the elastic constants of Vur- 
gaftman et alW^ and Ferhat and Zaouf^, respectively. It 
is worth noting that in the GaBia-Asi-a; case relaxed Ga- 
Bi bond lengths from our model follow the same trend re- 
vealed by x-ray absorption spectroscopy measurements^, 
where we see a decreasing Ga-Bi bond length in the pres- 
ence of Bi pairing and clustering, indicating increased 
localized strain in the alloy. 



TABLE I: Lattice constants, orbital self-energies, inter-atomic 
interaction parameters, spin-orbit couplings and valence band 
offsets for GaP, GaAs and GaBi in the sp 3 s* Hamiltonian. 
Lattice constants are given here in A; all other values are in 
eV. 





GaP 


GaAs 


GaBi 


a 


5.4505 


5.6532 


6.328 


E s>a 


-8.1124 


-8.6336 


-8.3774 


E s>c 


-2.1976 


-2.9474 


-5.6126 


Ep^a 


1.0952 


0.9252 


-0.1256 


Ep,c 


4.0851 


3.5532 


1.694 


E s *,a 


8.4796 


7.0914 


6.1262 


R * 

-C^s* ,c 


7.1563 


6.2 


5.8164 


ssa 


-1.8727 


-1.6835 


-1.3425 


ppa 


3.0989 


2.9500 


2.0003 


ppTY 


-0.7424 


-0.7420 


-0.6354 


SaP C (T 


1.8500 


2.3920 


2.3567 


ScPaU 


2.7312 


2.4200 


1.2025 


s* a p c a 


1.9998 


2.0400 


1.8051 


s*p a a 


2.1882 


1.7700 


0.5100 


\a 


0.0666 


0.405 


2.0142 


\c 


0.1734 


0.165 


0.1152 


VBO 


-0.5 
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TABLE II: 


Dimensionless scaling parameters, rj, for the sp 3 s* 


inter-atomic interaction parameters in GaP, GaAs and GaBi. 




GaP 


GaAs 


GaBi 


^tssa 


3.512 


3.512 


3.66 


T]sp<7 


4.3 


4.3 


4.08 


1]ppa 


3.2042 


3.2042 


2.2 


1Jpp-K 


4.2357 


4.2357 


3.24 


T]s* per 


5.5 


5.5 


5.5 


B. Tight-binding implementation 


of the band 




anti-crossing model 




The BAC model explains the extreme band gap bow- 


ing observed in Ga^^N^As 


i- x in terms of an interaction 


between two levels, one at energy Em associated with the 


extended CBE state \ip c ) 


of the Ga(In)As matrix, and 


the other at energy Ed associated with the localized N 


impurity defect states \ipjf 


), with the two states linked 


TABLE III 


VFF parameters 


and coefficient 


s for anharmonic 


corrections to the Keating potential for unstrained, bulk GaP, 


GaAs and GaBi. The VFF parameters ao and /3o are in units 


of N m . 


The anharmonicity constants A, B and C are 


dimensionless. 








GaP 


GaAs 


GaBi 


ao 


44.5 


41.49 


34.0 


Po 


10.69 


8.94 


5.0 


A 


7.2 


7.2 


8.2 


B 


7.62 


7.62 


7.62 


C 


6.4 


6.4 


6.4 
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by a N-composition dependent matrix element Vmd de- 
scribing the interaction between thenM The CBE energy 
of Ga(In)N ;c Asi_ :c is then given by the lower eigenvalue, 
E- , of the 2-band Hamiltonian 

/ Ed Vdm \ 

\V DM E m J 1 > 

A resonant feature associated with the upper eigen- 
value, E + , of Eq . (1) has been observed in PR 
measurements^2S12U j appearing in GaN^Asi-^ for x > 
0.2% and remaining a relatively sharp feature until x w 
3%, beyond which composition it broadens and weak- 
ens, when the resonant state becomes degenerate with 
the large -L-related conduction band density of states^. 

To investigate the resonant state \iPn), and its be- 
havior, we used the sp 3 s* TB Hamiltonian to cal- 
culate t he ele ctronic structure of ordered GaN^Asi-^ 
supercellsPM. By comparing the calculated host and 
alloy CBE states, \ip c ,o) and IV^i), in large supercells 
(Ga864NiAs863 and Gas64Ass64, respectively), we derived 
the nitrogen resonant state, \iPn), associated with an iso- 
lated N atom. In the BAC model, \ip c ,i) is a linear com- 
bination of IV'c.o) and |t/>at), with \iPn) then given by: 



As a result of the HH/LH degeneracy, the calculated al- 
loy VBE states each include an arbitrary admixture of 
the host matrix VBE states, IV^o)- We can therefore 
use the calculated alloy and host matrix TB wave func- 
tions to determine the four Bi resonant states \ipBi,i) as a 
linear combination of the GaBi^Yi-^ VBE, \ip v ,i,i)i and 
IVVcm)j the part of the host (GaY) matrix VBE with 
which the ith Bi state, \ipBi,i) interacts: 

4 

IVvm) - J2 \">Pi,o)(ipifi\ipv,i,i) 




in which the index I runs over the four-fold degenerate 
levels of the GaY VBE. 

In general, for an isolated Bi impurity, \4>Bi,i) exhibits 
Trf symmetry so it will be four-fold degenerate, with en- 
ergy FiBi- We can calculate explicitly Est and the in- 
teraction between IV'bm) and the GaY VBE, j3, by cal- 
culating the appropriate matrix elements of the full TB 
Hamiltonian for the Bi-containing alloy, H, for any one 
of the four degenerate Bi states: 



IV>, 



c,l) ~ 



O#c,o) 



(2) 



where, a — (V'e.olV'c.i) • The evolution of the conduction 
band structure in an ordered GaN^Asi-a, supercell can 
then be very well described by associating a localized 
resonant state with each N atom in the supercelP^l. 

In the dilute bismide case the BAC interaction is ex- 
pected to occur in the valence band. Previous analysis^ 
has included a BAC interaction with each of the three 
doubly degenerate valence bands, namely the heavy-hole 
(HH), light-hole (LH) and split-off (SO) bands. We find 
below that there is a BAC interaction with the HH and 
LH bands, but that there is no calculated BAC inter- 
action involving the SO band. For dilute bismides, the 
Hamiltonian analogous to (1) is therefore: 



Erh Vsi 
E LH 

V Bi E Bl 
V Bi 




(3) 



where E HH / LH are the (degenerate) host matrix VBE 
energies, E B i is the energy of the Bi-related resonant de- 
fect state and Vsi = P^/x is the composition-dependent 
interaction matrix element linking E B b/lh and Ebi- For 
E B i above (below) the host matrix VBE the interaction 
pushes down (up) the VBE. 

Because of the 4-fold degeneracy of the VBE states, 
we need to modify the approach used in Eq. (2) to de- 
termine the BAC interaction in GaBLjYi-^ (Y = P, As). 



E Bi = (ipBi\H\ip B i) 

V Bl (ipBi\H\t/) V fl) 







(5) 
(6) 



where \ip v ,o,i) = J2 l^z,o)(^.o|Vvi.i); with the sum again 

i=i 

running over the four-fold degenerate levels of the GaY 
VBE. 



C. Calculation of fractional T character spectra in 
alloy supercells 

In the ideal BAC model of Eqs. (1) and (3), the host 
matrix band edge states are expected to mix in the alloy 
with the localized Bi-related defect states. In order to 
explore this BAC interaction, we calculate Gr(E), which 
is the projection of the host matrix band edge states onto 
the full spectrum of levels in the alloy supercell. In what 
follows, we use the subscripts I, and fc, 1 to index the 
unperturbed host and Bi-containing alloy states, respec- 
tively. The Gr{E) spectrum is calculated by projecting 
the unperturbed G&n^n (Y = P, As) T states, \ipi,o), 
onto the spectrum for the GaArBijv/Yjv_Af alloy super- 
cell, {E k ,\i/j kA )}: 

g(Ei) 

Gr(£)^^/r, H T(E fc -£) (7) 

k 1=1 

fv,ki = |<V*,iM,o>| 2 (8) 



G 



Here we choose T (Ek — E) as a 'top hat' function of 
width 2 meV and unit height centered on Ek, and g(E[) is 
the degeneracy of the host (GaY) band having energy E\ 
at r (i.e. g{Ei) = 2,4 and 2 for the conduction, valence 
and SO bands, respectively). 

We use Gt(E) firstly to confirm the presence of a 
BAC interaction in ordered GaBi^-Pi^,,, and GaBi^Asi-^ 
structures, and then to analyse the effects of Bi pair and 
cluster states, as well as the evolution of the BAC in- 
teraction in disordered structures. We also use Gr(E) 
to explore the evolution of the conduction and SO band 
edges in GaBi^Asi-a; with Bi composition, showing that 
the variation in both cases is attributable to the conven- 
tional alloying effect. These calculations are discussed in 
detail in the next section. 



III. RESULTS AND DISCUSSIONS 

A. Bi state in ordered GaBi^Pi-^ and GaBLAsi 
supercells 



1. Gajv5iiPjv^i supercells 
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We have found previously that it is straightforward to 
identify and analyse a band anti-crossing interaction us- 
ing Eqs. (4) to (7) in ordered supercell structures where 
the localized defect state is in the energy gap, as was the 
case for GaN x Pi_j2H. It has been shown experimentally 
that an isolated Bi atom acts as an impurity in GaP, giv- 
ing rise to a defect level ~ 100 meV above the VBE 24 . 
Our tight binding model has been established to give a 
Bi defect state in the band gap of GaP, consistent with 
this experimental study. We therefore start our calcula- 
tions by investigating the electronic structure of ordered 
GajvBiiPjv-i supercells, in which a single substitutional 
Bi atom is inserted into the GaP matrix. Having es- 
tablished the main factors influencing the valence band 
structure in GaBi^Pi-a;, we then consider GaBi-rAsi-^, 
using the phosphide results to guide and support the in- 
terpretation of the arsenide calculations. 

Figure [3] plots the fractional T character, Gr{E), asso- 
ciated with the GaP VBE states in GajvBiiPAr_! super- 
cells as the supercell size, 2N (x = N" 1 ), decreases (in- 
creases) from (b) N = 2048 (x = 0.049%) to (f) N = 108 
(x = 0.926%). For reference, we also include in Fig. [3] (a) 
the calculated Gr(E) for a pure GaP supercell, in which 
100% of the HH/LH T character resides on the VBE. 

When we replace a single P atom by Bi in a 
Ga2048P2048 supercell, we see in Fig. [3] (b) that this gives 
rise to a four-fold degenerate defect level which lies above, 
and interacts with, the GaP VBE. The amplitude of the 
unperturbed GaP VBE states is distributed over two lev- 
els, corresponding to £V and E + , as expected from the 
band anti-crossing model of Eq. (3). We calculate us- 
ing Eqs. (4) and (5) that the isolated Bi impurity level 
lies 122 meV above the GaP VBE, and that each impu- 
rity state is highly localized, with 31% of its probability 



FIG. 3: Combined fractional HH/LH F character for a series 
of ordered GajvBiiPjv_i supercells. The zero of energy is 
taken at the GaP VBE. 



density localized on the Bi atom and the four nearest- 
neighbor Ga sites. We also calculate that the band anti- 
crossing interaction Vb% between the Bi-related defect 
levels, Ebi, and the VBE is of magnitude Vb% = f3\/x, 
with $ = 1.41 eV. 

The calculated magnitude of the interaction parameter 
/3 in GaBi^Pi-z is then comparable to that which we 
have calculated previously^ for GaN^Pi-^ (f3 = 1.74 
eV) and for GaN^Asi-^ (/3 = 2.00 eV). 

Having identified in this section that an anti-crossing 
interaction does indeed occur between Bi-related defect 
levels and the GaP VBE in GaBi a .P 1 _ a; , we now use this 
understanding to inform our analysis of the valence band 
structure of GaBi^Asi-a,. 



2. GajvBii Asn-i supercells 

Figure |4] plots the fractional T character, Gr{E), as- 
sociated with the GaAs VBE states in GajvBiiAs,/v_i 
supercells as the supercell size, 2N (Bi composition, x), 
decreases (increases) from (b) N = 2048 (x = 0.049%) 
to (f) N — 108 (x = 0.926%). For reference, we again in- 
clude as Fig. [4] (a) the calculated Gr (E) for a pure GaAs 
supercell, in which 100% of the HH/LH T character re- 
sides on the VBE. 

We first consider a single substitutional Bi atom in a 
4096 atom supercell (x = 0.049%, panel (b) of Fig. g). 
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FIG. 5: Plot of Bi resonant state energy E_b; and band anti- 
crossing parameter j3 as a function of the size of GaBi^Yi-a, 
(Y=P, As) supercell. Each supercell contains one Bi atom. 
The horizontal dotted line shows valence band edge energy in 
GaAs for GaBi^Asi-^ and in GaP for GaBi^Pi-^. 



We see for this case that the HH and LH Y character re- 
mains predominantly on the highest alloy valence states, 
with 98.2% of the Y character on these states. Of the 
remaining 1.8%, 1.3% is found on a single energy state 
w 117 meV below the GaAs VBE, with the other 0.5% 
spread over a series of lower lying levels. Using Eq. (4) , 
we identify that the state which mixes with each of the 
GaAs valence band edge states is a localized Bi resonant 
defect level, with very similar character to the localized 
defect state found in GaBi x Pi_a;. The Bi state is calcu- 
lated to be highly localized, with 36% of its probability 
density localized on the Bi atom and the four nearest- 
neighbor Ga sites. 

We calculate by evaluating Vm — (ipBi\H\ipv,o) for an 
isolated Bi impurity in a 4096 atom supercell that the 
band anti-crossing interaction Vbi between the Bi-related 
resonant defect level, Esi, and the GaAs VBE, E Vi g, is of 
magnitude Vb% = (3\/x, with (3 — 1.13 eV. The calculated 
magnitude of the interaction parameter j3 in GaBi^Asi-a, 
is then smaller than the value in GaBi 2 ,P 1 _ :I ., reflecting 
that there is less difference in electronegativity and size 
between a Bi and an As atom than there is between a Bi 
and a P atom. 

When the size of the supercell is decreased from 4096 
atoms to 216 atoms (panels (b) to (f) in Fig. [4|, the 
interaction of the Bi-related resonant defect state with 
the GaAs VBE increases, thus pushing the valence band 
edge, E_, upward in energy. The lower BAC level is 
pushed downward in energy and its fractional Y charac- 
ter, Gr(E), also increases as a result of the increasing 
strength of the interaction. Overall, the Y character as- 
sociated with the lower BAC level is in all cases spread 
over several supercell levels (clearly visible in panels (b), 
(c) and (d), but running out of scale to lower energies 
in panels (e) and (f)), reflecting that the Bi states are 
resonant with a high density of host matrix valence band 
states. 

Figure [5] shows the calculated values of f3 and Eb% 
for a series of ordered 2N atom GaArBiiYAr_i (Y = P, 
As) supercells, where N is increased from 32 to 4000 
atoms, retaining a cubic supercell at each N. Both the 
GaBi^Asi-a, and GaBi^Pi-a, supercells exhibit similar 
trends. As the supercell size increases, the calculated 
strength of the band anti-crossing interaction between 
the Bi impurity state and the GaAs VBE reduces. For 
larger supercells (N > 2000 atoms) we approach the di- 
lute doping limit and the trend stabilises. We also find 
that the calculated energy of the Bi resonant states, Ebu 
increases with increasing supercell size. The Bi impurity 
state is calculated to lie 607 meV below the GaAs VBE in 
a Ga32BiiAs3i supercell, and its value increases to —0.18 
eV in a Ga4oooBiiAs3gg9 supercell. This change arises 
due to the evolution of the GaY (Y = P, As) supercell 
zone center band structure with increasing supercell size. 
The GaY supercell zone center states form a complete 
set of basis states from which to construct the localized 
Bi states^. For the smallest (64 atom) supercell, the 
highest valence levels folded back to the zone center lie 



(a) Bi pair (b) Bi cluster type 1 (c) Bi cluster type 2 




FIG. 6: Three different Bi nearest-neighbor environments are used in our study of pairing and clustering effects: (a) a single 
Bi 'pair' (i.e. an isolated Ga atom having two Bi nearest neighbors), (b) a cluster of Bi atoms is placed inside the supercell 
such that a single Ga atom has three Bi nearest neighbors ('cluster type 1'), and (c) a cluster of Bi atoms containing three Bi 
atoms in the form of a Bi-Ga-Bi-Ga-Bi chain. In this configuration, two Ga atoms have two Bi nearest neighbors ('cluster type 
2'). Note that the nearest-neighbor environments are depicted here in their unrelaxed configurations. The bonds in the red 
color indicate the Bi pair and cluster connecting bonds. 



1.45 eV (1.30 eV) below the VBE in GaAs (GaP). As the 
supercell size increases, more states get folded back with 
energies closer to the VBE. These states can be used to 
construct the Bi localized state energy, thus accounting 
for the calculated upward shift in the energy Eb% with 
increasing supercell size. Based on the trends shown in 
Fig. [5j we predict values of approximately —0.18 eV and 
1.13 eV for Est and f3 for GaBi 2 .As 1 _ :! . alloys in the dilute 
doping limit, respectively. 



B. Bi pair and cluster states 

Recent experimental investigations using x-ray absorp- 
tion spectroscopy^ indicate that Bi atoms in dilute 
GaBi^Asi-z alloys tend to form localized pairs and clus- 
ters, even at relatively low Bi compositions. The pres- 
ence of such pairs and clusters of N atoms has been ob- 
served in earlier studies on GaN x Asi_ x alloys with dilute 
N concentrations^ and it has been shown that such pairs 
and clus ters s ignificantly affect the electronic structure of 
the allo} * 20 * 46 4 We investigate here the impact of Bi pairs 
and clusters on the electronic structure of GaBi^Asi^a;. 

We compare the electronic structure of seven dif- 
ferent nearest-neighbor environments in a 4096 atom 
supercell: (i) a Ga 2 o48As 2 o48 supercell (Bi-free), (ii) 
Ga2048BiiAs2047 containing an isolated Bi atom, (hi) 
a Ga2048Bi2As2046 supercell containing two isolated Bi 
atoms (where each Bi atom has no Bi belonging to ei- 
ther its first, second or third shell of neighboring atoms), 

(iv) a Ga2048Bi2As2046 supercell containing a pair of Bi 
atoms (i.e. a Ga atom having two Bi nearest neighbors), 

(v) a Ga2048BiaAs2045 supercell containing three isolated 
Bi atoms (with each Bi atom satisfying the same 'isola- 
tion' criterion as in (ii)), (vi) a Ga204sBi3As2045 supercell 



containing three Bi atoms forming a cluster, such that a 
single Ga atom has three Bi nearest neighbors (we shall 
henceforth refer to this nearest-neighbor environment as 
a cluster of type 1), and (vii) a Ga204sBi3As2045 super- 
cell containing three Bi atoms forming a Bi-Ga-Bi-Ga-Bi 
chain, with the two Ga atoms in the chain each having 
two Bi nearest neighbors (this will be referred to as a 
cluster of type 2). The nearest-neighbor environments 
(iv), (vi) and (vii) are depicted in Figs. [6] (a), (b) and 
(c), respectively. 

From the relaxation of the supercell configurations de- 
scribed in (ii) to (vii), we calculate the average relaxed 
Ga-Bi bond length in each case. For a single Bi atom, 
the relaxed Ga-Bi bond length is 2.647 A. When we intro- 
duce a Bi pair, the Bi-Ga bond length to the shared Ga 
atom reduces to 2.640 A. A further reduction in the bond 
length occurs for Bi clusters. For example, the average 
Bi-Ga bond length to the shared Ga atoms is 2.640 Afor 
a cluster of type 2. These trends are similar to the ones 
reported by experimental studjEl, where a reduction in 
the bond lengths is observed for pairs and clusters. The 
x-ray absorption measurements indicated average bond 
length d av values of 2.627 A, 2.621 A, and 2.619 Afor 
isolated, pairs, and clusters of Bi atoms. 

Figure [7] shows the calculated fractional GaAs VBE T 
character for all seven of these nearest-neighbor environ- 
ments. We recall in panel (b) that inclusion of a single 
Bi atom results in the mixing of a small amount (« 2%) 
of the GaAs VBE T character into alloy states in the 
VB. Similar weak interactions can be seen in panels (c) 
and (e), for the cases where two and three isolated Bi 
atoms are incorporated in the alloy. In these cases, a 
large fraction (96% in panel (c) and 95% in panel (e)) of 
Gr(E) again resides on the highest alloy valence state, 
with the remaining 4-5% distributed over a number of 
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lower valence levels. 

Turning to the Bi pair in panel (d), the T character 
is in this case distributed over four states, which form 
two interacting, anti-crossing pairs (indicated by red and 
black bars) . One of the pairs shown with the black color 
bars (w 5.3% weight at -90 meV and w 94% at 10 meV) 
behaves in a very similar way to the pair of states asso- 
ciated with isolated Bi atoms in Fig. 0(b) and (c), while 
the second pair (red bars) behaves quite differently, with 
i=3 55% r character on the lower state at -27 meV and f=s 
44% T character on the upper state at 43 meV, associ- 
ated with a Bi pair-related defect level just above the host 
matrix VBE. We can understand the behavior of these 
defect states based on interactions between the resonant 
states associated with two isolated Bi atoms. The forma- 
tion of the Bi pair breaks the Td symmetry about each Bi 
site, lifting the four-fold degeneracy of the Bi-related lo- 
calized states to give a pair of two-fold degenerate states, 
one of which has a small interaction and the other a large 
interaction with Bi states on the neighboring Bi site. The 
interaction between the two Bi states with large overlap 
then gives rise to the defect level which we calculate to be 
in the g ap, s imilar to the behavior of the N-N pair state 
in GaAgSEU, while the interaction between the other two 
sets of Bi states gives the resonant level calculated to be 
90 meV below the host matrix VBE. 

A similar analysis holds for the two types of Bi cluster 
considered. The type 2 cluster (chain of 3 Bi atoms) has 
a very similar spectrum to the Bi pair, with a stronger 
interaction between the lower defect and the VBE states 
than in the pair case, and with the higher defect state 
being slightly deeper in the energy gap. We estimate 
that the energy level of the Bi defect state (Esi) is 10 
meV, 8 meV and 17 meV above the GaAs VBE for the Bi 
pair and the clusters of type 1 and 2, respectively. The 
lower defect state is shifted to -48 meV in the presence of 
a type 1 cluster, where one Ga atom has 3 Bi neighbors. 
Overall, the results follow similar trends to those found 
previously for N-related defect clusters in GaAs, where 
the defect levels tended to move deeper into the GaAs 
energy gap with increasing cluster size. 

Table IIVI summarises the calculated effect of the dif- 
ferent Bi configurations on the lowest conduction band 
energy, E Cl the top two valence bands, E v \ and E V 2, and 
the spin-split-off band, Ego- Also shown in parentheses 
below each energy are the calculated change in energy for 
the band (per % Bi incorporation, relative to the corre- 
sponding GaAs band edge) as well as the overlap of the 
state with the corresponding unperturbed GaAs state. 
We see from these calculations that E c and Eso ar e both 
largely insensitive to the nearest-neighbor environment 
in the Bi-containing supercells. The conduction band ex- 
hibits a decrease of 28 meV per % Bi, which is linear in Bi 
composition, and retains a 99.9% overlap with the GaAs 
CBE for all Bi configurations considered. Thus, the vari- 
ation in the conduction band and spin-split-off energies 
with increasing Bi composition can be well explained in 
terms of conventional alloy models, irrespective of disor- 
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FIG. 7: Combined fractional HH/LH V character for a series 
of disordered Ga204sBiAf AS204S-M supercells, for several dif- 
ferent nearest-neighbor environments. The zero of energy is 
taken at the GaAs VBE. 

der in the form of Bi pairs and clusters in the alloy. We 
note that unexpectedly large values have been deduced 
for the conduction band edge mass in GaBi^Asi-j; alloys 
for x < 6%^. Further studies are required to identify 
the cause of the anomalous CBE mass values deduced 
experimentally: our calculations indicate that it should 
be possible to describe the CBE mass in randomly disor- 
dered GaBi^Asi-a, using conventional alloy models. 

C. Electronic structure of disordered 

Ga 2 048Bi M P 2 048-M supercells 

Having established the character of isolated Bi single, 
pair and cluster impurity states in sections III-A and III- 
B, we now investigate the electronic structure of ran- 
domly disordered GaBi^Pi-a; alloys, as we increase the 
number (M) of Bi atoms in a Ga204sBijif P2048-M super- 
cell. In order to trace the evolution of the highest valence 
states, we diagonalise the sp 3 s* Hamiltonian for the su- 
percell and use the resulting alloy states to calculate the 
fractional V character, G T (E), for the VBE (Eqs. (7) 
and (8)). 

The calculated Gr(E) spectra for a range of Bi com- 
positions are shown in Fig. [81 We see in panel (b) for 
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FIG. 8: Combined fractional HH/LH V character for a se- 
ries of disordered Ga204sBiAfP2048-A/ supercells, for several 
Bi compositions, x = ^ffg • The zero of energy is taken at the 
GaP VBE. 



a single Bi atom in the supercell (M = 1, x = 0.05%) 
that the amplitude of the unperturbed GaP VBE states 
is distributed over two levels, corresponding to £L and 
E+, as previously discussed in section III-A. 

In the supercells used in calculations containing M > 
10 Bi atoms (x > 0.49%) we see that the calculated split- 
ting between the and E + levels does not increase 
smoothly with increasing Bi composition. In addition, 
because of the reduced symmetry in these supercells and 
the inclusion of Bi pairs and clusters, the two highest 
valence states are no longer degenerate and we see an 
overall broadening of the spectra. 

The pair and cluster states introduce higher energy 
defect states, lying above the isolated Bi states in the 
energy gap, as can be seen at higher Bi composition in 
Fig. [8] The Bi defect state distribution in CaBi^Pi^^ 
then has a similar form to that observed in CaN^Pi-^, 
where isolated N defect states are just below the GaP 
CBE, with defects states associated with N pairs and 
other cluster levels lying deeper in the energy gapP^. 

We can see from Fig. [8] that there is a general anti- 



crossing interaction between the valence band states and 
the Bi defect states in the gap, with an increasing sepa- 
ration between the VB-related states (below eV) and 
the Bi states in the gap (> 0.2 eV) with increasing Bi 
composition. However, because the gap defect levels are 
distributed over a range of energies, it is not possible to 
describe the band structure close to the energy gap using 
the two-level band anti-crossing model of Eqs. (1) and 
(3). We showed previously for GaN^Pi-^ that it is possi- 
ble to treat the anti-crossing interaction between the GaP 
host matrix T CBE and a distribution of N-related de- 
fect levels by treating explicitly the interaction between 
the GaP T conduction state and the full distribution of 
N defect levelfP. We see from Fig. [8] that this should 
also be the case in GaBi 2 .P 1 _ :I .. When the supercell only 
included isolated Bi atoms then a clear anti-crossing in- 
teraction was observed in Fig. [3j in agreement with Eq. 
(3). However, when the calculation includes a wider dis- 
tribution of Bi environments (which is expected experi- 
mentally even for x < 1%) then the 2-level BAC model 
is no longer able to describe the observed evolution of 
the VBE states. We expect that a model similar to that 
used for GaN^Pi-^ could also be used to describe the 
evolution of the band structure of GaBi^Pi-^. We do 
not however pursue that further here. 



D. Electronic structure of disordered 

Ga 2( )48Bi A /As2048-M supercells 

1. Evolution of GaBi x As±- x valence states with x 

Figure [9] shows the calculated evolution of 
Gr(E) for the highest valence states in a series of 
Ga2048Bij\/As2048-A/ supercells, where the supercell 
energy spectrum is in each case projected onto the un- 
perturbed GaAs valence band edge heavy-hole (HH) and 
light-hole (LH) wave functions (see Eqs. (7) and (8)). 
To enable comparison with the results of our GaBi^Pi-a; 
calculations, we consider identical Bi distributions 
within the supercells, with all P atoms replaced by As. 
For reference, the lowest panel, (a), shows the calculated 
Gt{E) for an unperturbed Ga204sAs2048 supercell, which 
is simply a single peak of unit amplitude, centered at 
the location of the Ga 2 o4sAs2048 HH/LH energy levels. 

We observe in Fig. [9] that the calculated VBE energy 
shifts upwards with increasing Bi composition, x, and 
that the overall T character of the valence band maximum 
state decreases with increasing x, as expected from the 
BAC model of Eq. (3). The solid circles and triangles in 
Fig. ~ 

of each of the two highest valence 



10 show how the calculated fractional V character, 

fr.v = KVvolVvi)! 2 , 
states decreases with increasing Bi composition in the 
supercells considered. 

The solid line in Fig. [10] shows the expected variation 
in T using Eq. (3) with: E HH / LH = 1.03 x eV, E B i — 
-0.18 eV and V Bi = /3,/x, with = 1.13 eV. It can be 
seen that a good fit is obtained at low x, with the quality 
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TABLE IV: Calculated effect of Bi isolated, pair, and clusters of atoms on the lowest conduction band, the highest two valence 
bands, and the spin-split-off band edges of Ga204sBi.M AS2048-M supercells, including energy of given state and in parentheses 
the calculated change in energy for the band (per % Bi incorporation, relative to the corresponding GaAs band edge) as well 
as the overlap of the state with the corresponding unperturbed GaAs state. 



Alloy Composition 


x(%) 


Ed (eV) 
(meV/%Bi, %) 


E„i (eV) 
(meV/%Bi, %) 


E„ 2 (eV) 
(meV/%Bi, %) 


E so (eV) 
(meV/%Bi, %) 


Unperturbed GaAs 





1.519 

(0, 100) 


0.000773 

(0, 100) 


0.000773 

(0, 100) 


-0.3518 

(0, 100) 


Single Bi-atom 


0.049 


1.5176 

(28.57, 99.97) 


0.00456 

(77.28, 98.44) 


0.004359 

(73.18, 98.29) 


-0.3521 

(7.06, 99.95) 


2-isolated Bi-atoms 


0.098 


1.5163 

(27 55 99 93") 


0.008106 

(74 82 96 92*1 


0.007746 

(71 15 97 15^1 


-0.35244 

(6 53 99 94~1 


Bi-pair 


0.098 


1.5162 

(28.57, 99.93) 


0.04298 

(430.68, 44) 


0.009488 

(88.92, 94.02) 


-0.35244 

(6.62, 99.34) 


3-isolated Bi-atoms 


0.146 


1.5149 

(28.08, 99.9) 


0.01146 

(73.2, 95.8) 


0.01078 

(68.51, 96.44) 


-0.35273 

(6.33, 99.75) 


Bi-cluster Type 1 


0.146 


1.5148 

(28.76, 99.88) 


0.09803 

(666.14, 18.14) 


0.0223 

(147.44, 77.73) 


-0.3526 

(5.6, 99.75) 


Bi-cluster Type 2 


0.146 


1.5149 

(28.08, 99.9) 


0.05117 

(345.18, 41.06) 


0.01493 

(96.96, 89.34) 


-0.35273 

(6.35, 99.64) 



of the fit decreasing at higher values of x, as we move 
into a disordered regime in the alloy. 

A resonant feature associated with the upper eigen- 
value, E +1 of the Hamil tonian (1) has been observed 
in PR measurement d 26 ' 27 1 on GaN^Asi-^ for x > 0.2%, 
giving a relatively sharp feature until x reaches approxi- 
mately 3%, beyond which composition the feature broad- 
ens and weakens, when the resonant state becomes degen- 
erate with the L-related CB levels^. It has not proved 
possible to observe such a feature at any composition 
in GaBi^Asi-a;, and this has led to controversy as to 
whether there is a VB-BAC interaction present in the 
allojfl 3 . It can be seen from Fig. [9] that the GaAs VBE 
character in GaBi^Asi-a; is mixed into a wide range of 
the alloy valence states so that, unlike in the GaBi a; Pi_ x 
case, no single state can be identified as the E + level in 
GaBi^Asi-a;. This strong mixing is to be expected, be- 
cause the width of t he re sonant level is proportional to 
the density of states^EH. The large density of valence 
band states then explains the failure to observe an E + 
level in GaBi^Asi-^ PR spectra. 



2. Band gap and spin- orbit- splitting energies versus x 

Figure [TT] shows how the fractional T character of the 
SO band edge evolves with increasing Bi composition. 
As the Bi composition increases, the T character spreads 



out over several states and it becomes unclear to choose 
a particular energy for E SO - The values of the E so en- 
ergies from panels (b) to (g) of Fig. 11 are obtained by 
taking a weighted average of the energies with their cor- 
responding r character as well as by picking the energies 
with the highest T character. The values from the two 
methods come out to be very close (within 0.1 meV of 
each other) and we deduce that both methods are equally 
good in this case. The Bi-induced disorder leads to a 
broadening of the projected T character, Gr(E), with in- 
creasing Bi composition, but no anti-crossing behavior 
is observed in this broadened spectrum. Hence, the shift 
and broadening of the SO state projection is attributable 
to conventional alloying and disordering effects. For this 
reason we conclude that our calculations do not support 
the existence of a Bi resonant state interacting with the 
SO band, as introduced in the kp Hamiltonian of Alberi 
et al¥*. 

Figure [12] shows the calculated variation of the lowest 
conduction, the two highest valence, and the spin-split- 
off band edge energies with increasing Bi composition x 
in the 4096 atom GaBi x Asi_ x supercells. It can be seen 
that the calculated conduction band edge energy shifts 
down linearly by about 28 meV per % of Bi replacing As. 
The solid squares in Fig. 10 show that E c retains > 97% 
GaAs conduction band T character for x < 5%. With a 
linear decrease in the energy of E c and the T character 
remaining close to 100%, we conclude that the evolution 
of the conduction band states is also best described using 
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FIG. 9: Combined fractional HH/LH V character for a series 
of disordered Ga2048BiwAs2048-M supercells for several Bi 



compositions, x 
GaAs VBE. 



The zero of energy is taken at the 



a conventional alloy model. 

An important feature of GaBi^Asi-^ alloys is that, as 
the Bi composition increases, the decreasing band gap 
Eg = E c - E v i and the increasing spin-orbit splitting en- 
ergy Ago = E v i — Eso results in a crossing of these 
two energies, with E g < Aso a t higher Bi composi- 
tions. It has been suggested that this should lead to 
a suppression of the dominant CHSH Auger recombi- 
nation mechanism in GaBi^Asi^^ at higher Bi compo- 
sitions, which could be highly beneficial for the design 
of high efficiency lasers and semiconductor optical am- 
plifiers (SOAs) operating at telecommunication and at 
longer wavelengths^. Recent PR measurements^ have 
confirmed a crossing of the energy gap and spin-orbit- 
splitting energies in bulk GaBi^Asi-a; alloys for compo- 
sitions above about £=10.5%. 



Figure 13 plots the calculated band gap E g and spin- 



orbit splitting Aso as a function of Bi composition in 
4096 atom GaBi-rAsi-^ supercells. The calculations were 
undertaken using the low temperature GaBi^Asi^j, band 
parameters presented here (£' ff (GaAs)= 1.519 eV), and 
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FIG. 10: Plot of fractional F character, /r (see Eq. 8), of the 
lowest conduction band (E c ) and the highest two valence band 
(E„i and E„2) edges as a function of Bi composition. The 
solid line shows the valence band fr character calculated using 
parameters in Eq. 3 extracted from tight-binding calculations 
of supercells containing isolated Bi atoms. 
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FIG. 11: Combined fractional SO T character for a series 
of disordered Ga204sBiMAs2048-M supercells for several Bi 
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FIG. 12: The lowest conduction band edge E c , the highest 
two valence band edges E v i, E V 2, and the spin-split-off band 
E so energies are plotted as a function of Bi composition in 
GaBi x Asi-rc supercells. The conduction band E c edge de- 
creases with a linear slope, changing by about 28 meV for 
every 1% increase in the Bi composition. The highest two va- 
lence band edges E v i and E V 2 exhibit a non-linear change. 
For small Bi composition, the slope is large 53 meV/1% 
Bi. For large Bi compositions, the slope decreases to about 
30 meV/1% Bi. The spin-split-off band E so also decreases 
following a nearly linear slope of 5.5 meV/1% Bi. 



the CBE energy has then been shifted down by 92 meV 
for all Bi compositions, to fit to the GaAs energy gap at 
room temperature. Each data point shown in the figure is 
calculated by taking the average of three different random 
configurations of GaBi-rAsi-^. The results presented in- 
dicate that the decreasing band gap E g and the increas- 
ing spin-orbit splitting do indeed lead to a crossover of E g 
and Aso at about 10.5% Bi composition, in agreement 
with PR measurement^^. For a broader comparison with 
other experiments, we have also included four sets of ex- 
perim entally m easur ed b and gap and spin-orbit splitting 
values^* 10 * 12 ! in Fig. 13 along with our calculated values. 
Our results exhibit a close match with the experimental 
data points over the entire composition range considered, 
confirming the validity of our model. 
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FIG. 13: The band gap E g and spin-orbit-splitting Aso en- 
ergies plotted as a function of Bi composition in 4096 atom 
GaBizAsi-a, alloy supercells. The solid circles and squares 
show the calculated data from our model. The triangles and 
open circles indicate the measured PR data from experiment 
(Z. Batool et aZ.^). Four other sets of experimental data are 
also included for broader comparison. 



IV. CONCLUSIONS 

In summary, we have developed a model based on the 
valence force field (VFF) and sp 3 s* tight binding meth- 
ods to study the dilute bismide alloys GaBi^Asi-^ and 
GaBi^Pi-a;. Our results confirm that the evolution of the 
highest valence states in GaBi a; Yi_ 2 ; (Y = P, As) can be 
described using a band anti-crossing model. This has, 
until now, proved difficult to confirm experimentally be- 
cause one of the key BAC features, the E + state observed 
in GaN-rAsi-z, is absent in GaBi^Asi-^. We showed 
that this higher-lying E + level which has been identi- 
fied in GaNzAsi-z cannot be observed in GaBia;Asi_ x , 
due to the resonant levels being strongly broadened be- 
cause of their interaction with the large density of va- 
lence states. In addition, we find that the presence of Bi 
pairs and clusters in a randomly disordered alloy intro- 
duce further defect levels, which also contribute to the 
suppression of the E + feature. We calculate that the ex- 
perimentally observed strong band gap reduction is only 
partially explained through the valence band BAC inter- 
action. We find that the conduction band edge energy 
also decreases with increasing Bi composition, and that 
both band edge shifts contribute approximately equally 
(at large Bi compositions) to the reduction in the band 
gap energy of unstrained GaBi^Asi-j; alloys. The com- 
puted dependence of band gap on Bi composition closely 
matches with experimental measurements. We calculate 
that the band gap and spin-orbit-splitting energies cross 
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each other at 10.5% Bi composition in GaBi x Asi_a;, 
consistent with recent experimental measurements^. It 
has been suggested that this crossover may eliminate the 
dominant CHSH Auger recombination loss mechanism in 
longer wavelength photonic devices. Overall, we conclude 
that dilute bismide alloys open up new opportunities for 
band structure engineering and its application, with sig- 
nificant potential benefits for GaAs-based optical devices. 
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